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Abstract. 

A hydro dynamic formulation for accretion flow channeled by a dipo- 
lar magnetic field is constructed using a curvi-linear coordinate system 
natural to the field structure. We solve the hydrodynamic equations and 
determine the velocity, density and temperature profiles of the post-shock 
accretion flow. The results are applied to accretion flows in intermedi- 
ate polars. We have found that for systems with massive white dwarfs 
(~ 1 M0) the temperature profiles in the flow can differ significantly to 
those obtained from models in which the accretion column is assumed to 
be cylindrical. 



1. Introduction 

In magnetic cataclysmic variables (mCVs), the accretion flow near the white 
dwarf is channeled by the magnetic field. The flow is initially supersonic but 
becomes subsonic near the white-dwarf surface, and an accretion shock is formed. 
The shock heats the accreting material to temperatures of ~ lO'^ — 10^ K. The 
shock-heated material is ionized and emits bremsstrahlung X-rays and cyclotron 
optical/IR radiation (Lamb & Masters 1979). The height of the shock above the 
white-dwarf surface is determined by the cooling processes. It had been believed 
that the height is small 10^ — 10^ cm, see Wickramasinghe & Meggitt 1985), 
and therefore the shock-heated emission region could be treated as a stack of thin 
pancake-like plasmas on the white-dwarf surface. More recent studies showed 
that for some mCV parameters the shock height can be comparable to the 
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Figure 1. (Left) The {u,w) plane of the curvi-hnear coordinate sys- 
tem used in the formulation. Solid curves represent dipole magnetic 
field lines, and dashed curves represent field equi-potential surfaces. 
(Right) The w component of the gravity force, 5^, as a function of 
height h, for field lines with different foot points. The force component 
is normalized to that of the field line along the polar axis, g, and the 
height is normalized to the shock height H. The lines (from top to 
bottom) correspond to the gravity force along the field lines with foot 
points at magnetic co- latitude 9 = 0°, 6°, 12° and 18° respectively. 



white-dwarf radius when the effects of gravity along the flow are considered (see 
e.g. Cropper et al. 1999). The latter studies predicted X-ray spectra in good 
agreement with the grating spectra of mCVs obtained by XMM-Newton and 
Chandra (see e.g. Cropper et al. 2002). 

Most existing hydrodynamic calculations of the post-shock accretion flow 
in mCVs assume a simple geometry, such as a semi-infinite slab or a cylinder 
(Aizu 1973; Chavelier & Imamura 1982; Wu, Chanmugam & Shaviv 1994; Crop- 
per, Ramsay &: Wu 1998). These simple geometries are good approximations, 
provided that the shock height is small and the flow is azimuthal symmetric and 
aligned close to the polar axis. For accretion flows in intermediate polars, these 
two conditions are not always satisfied, and hence, the effect of the magnetic- 
field geometry needs to be investigated in order to obtain accurate velocity and 
temperature profiles of the accretion flow. 

Here, we present a hydrodynamic formulation in terms of a curvi-linear co- 
ordinate system natural to channeled flows in dipole magnetic fields. We solve 
the hydrodynamic equations and determine the velocity, density and tempera- 
ture structures of the post-shock region. The results are compared with those 
of Cropper et al. (1999), in which a azimuthal cylindrical accretion column is 
assumed. 



2. Formulation 



We assume the accreting material is an ideal gas, with a pressure P, density p 
and temperature T, related by P = pk^T/ fim^i, where is the Boltzmann 
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constant, toh is the hydrogen mass, and ii is the mean molecular weight. In 
stationary accretion, the hydrodynamic equations governing the flow are 

V-(pv) = 0, 
(v-V)v+^VP = g, 

(v.V)P-^(v.V)p = -(7-l)A, 

where v is the velocity of the flow, 7 is the adiabatic index of the gas, g is 
the gravitational acceleration, and A is the effective cooling function. (See Wu 
(2000) and references therein for the hydrodynamic formulation of accretion flow 
in mCVs.) 

We consider a coordinate system {u, w, cp), in which the orthogonal unit 
vectors ii, w and (p are defined as follows: u is on the equi-potential surface of 
the field at a fixed azimuthal angle, w is along the magnetic field line, and ip is 
the same as that for the azimuthal coordinate in the spherical coordinate system 
(see Fig. 1, left panel). The metrics of the coordinate system, hi, h2 and /13, are 
given by 



^3 = 1 hr- 




hi = 

ho = W(^l +(^1 +(^ 



where {x,y,z) is the Cartesian coordinate system. We may then express the 
hydrodynamic equations in terms of (u, w, ip). 

We integrate the mass continuity equation directly and obtain hih^pv = 
C, where C is a constant. Substituting this into the other two hydrodynamic 
equations yields 




Here, ^ = v + {P/pv), is the gravitational acceleration projected on w (Fig. 1, 
right panel), and the function H{u,w) = d\ii{hihz) /dw, describes the change of 
cross-section area of the flow-flux tube. 

Our cooling function A consists of two parts, one for bremsstrahlung cooling 
and another for cyclotron cooling. We adopt the composite cooling function of 
Wu (1994) and Wu et al. (1994) (see also Cropper et al. 1998). As cyclotron 
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Figure 2. The quantity pv as a function height along the field lines, 
normalized to the values, (pf)wd) at the white-dwarf surface. From top 
to bottom, the foot points of the field lines are at magnetic co-latitudes 
6* = 0°, 6°, 12° and 18° respectively. 

cooling is practically unimportant in intermediate polars, we simply set the 
efficiency parameter eg = 0, i.e. A = Ap^i^P/ p)^/"^ , where the constant A = 
3.9 X 10^^ in c.g.s. units. We assume a cold stationary wall for the lower 
boundary condition as in Chevalier & Imamura (1982) and Wu et al. (1994). 
For the upper boundary condition a strong adiabatic shock is considered, but 
with this coordinate system its location must be determined self-consistently. 
Detailed treatment of the upper boundary condition will be presented in Canalle 
et al. (2003). With the cooling function and the boundary conditions defined, we 
can now integrate the hydrodynamic equations along the field lines numerically 
and obtain the velocity, density and temperature profiles of the fiow. 

3. Results and Discussion 

In all calculations, unless otherwise stated, the white-dwarf mass Mwd = IMq 
and the specific mass accretion rate m = 2.0 g cm~^s~^. The white-dwarf radius, 
R^Ai is determined by the Nauenberg (1972) mass-radius relation. 

For comparison we also calculate the velocity, temperature and density pro- 
files of a flow with an azimuthal symmetric cylindrical accretion column using 
the formulation in Cropper et al. (1999). The two main differences in the two 
formulations are that for flows channeled strictly by dipolar magnetic flelds the 
cross section of the flow-flux tube decreases with distance to the white-dwarf 
surface and that the gravitational acceleration varies in strength and direction 
along the flow. As shown in Figure 1 (right panel), for the parameters assumed 
the gravitational accelerations are smaller for the flows along field lines with 
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Figure 3. The velocity profiles of the post-shock accretion flow. The 
solid, dashed, dash-dotted and dotted curves correspond to flows along 
fleld lines with foot points at magnetic co-latitudes 9 = 0°, 6°, 12° 
and 18° respectively The dash-dot-dot-dotted curves curve (labeled 
CWRK) correspond to the velocity proflle obtained from the calcula- 
tion following Cropper et al. (1999). The curves are terminated at the 
heights where the shocks are located. The velocities of the curves are 
shifted progressively upward each by 0.1 x 10*^ cm s^^ to amplify their 
differences. 

foot point co-latitudes 9 farther away from the pole. With the parameters as- 
sumed in Figure 1, the value of 5«;/|g| at the shock is about 0.985 for 9 = 18°: 
the value can be much smaller for larger 9 and also for parameters which yield 
shock heights signiflcantly above the white-dwarf surface. Figure 2 shows the 
pv proflles for various 9. For small 9, the proflles are similar, but when 9 is 
substantially large the deviation becomes signiflcant: compare the proflles of 
9 = 12° and 18° with those of smaller 9. 

Figure 3 shows the velocity proflles of the flows. The shape of the proflles 
for the flows are quite similar regardless of the differences between the values of 
9. However, the shock heights are different, and hence the shock temperatures 
are different accordingly. Interestingly, the proflles are quite different to those of 
the cylindrical accretion column (Cropper et al. 1999). In particular, the latter 
yield a smaller shock height, thus a hotter and more compact post-shock emission 
region. This will affect the properties of the optical/IR cyclotron radiation as 
they are emitted mainly from the hotter part of post-shock region. The effects 
on the bremsstrahlung X-ray continuum may be less affected as it is emitted 
from the cooler bottom of the post-shock region. 

Figure 4 shows the temperature proflles for systems with various white- 
dwarf masses. Both the flows in cylindrical columns and flows conflned by dipo- 
lar flelds show that the shock temperatures increases with the white-dwarf mass. 
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Figure 4. Temperature profiles of flows along a dipolar field in com- 
parison with flows assuming a cylindrical accretion column (labeled 
CWRK, Cropper et al. 1999) for white-dwarf masses 1.0, 0.9, 0.8 and 
0.7 Mq (curve pairs from top to bottom). For the flows along a dipolar 
field, the foot point of the field line is fixed at the magnetic co-latitude 
Q = 18°. In all cases, the mass accretion rate fn = 2.0 g cm~^s~^. 



For white-dwarf with masse ^0.7 Mq, the profiles for the two cases are similar, 
with the model with cylindrical columns giving slightly higher temperatures. 
However, when the white-dwarf mass is sufficiently large (~ 1.0 Mq), the differ- 
ences between the two cases becomes obvious. Also, the model with cylindrical 
columns tends to over predict the gas temperatures. The difference between the 
two cases can be explained qualitatively as follows. Larger white-dwarf masses 
result in higher shock temperatures and hence larger shock heights above the 
white-dwarf surface. When the shock height becomes too large, the cylindri- 
cal accretion column is not a very good approximation, and the variation in the 
gravitational acceleration along the flow becomes significant. The field geometry 
has stronger effects on the accretion onto massive white dwarfs than low-mass 
white dwarfs. Therefore, when modelling the X-rays and optical/IR radiations 
from the shock-heated accreting gas in the systems with massive white dwarfs, 
one must consider the effects due to the field geometry. 

These finding may explain the finding of Cropper et al (1999) and Ramsay 
(2000) who found that the mass of the white dwarf in mCVs were biased towards 
higher masses compared to isolated white dwarfs when they fitted X-ray spectra 
with the cylindrical formulation. Work is currently in progress to implement the 
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affect of the field geometry and determine how it will affect the resulting white 
dwarf masses. 

In summary, we present a hydrodynamic formulation using a curvi-linear 
coordinate system defined by the dipole field configuration. We use this formu- 
lation to calculate the velocity, density and temperature profiles of the accretion 
flows in mCVs. We have shown that effects of field geometry are important for 
mCVs with high-mass white dwarfs. Our formulation improves the calculations 
of temperature and density structures of post-shock emission region in mCVs, 
in particular, the intermediate polars, and thereby predict more accurately the 
X-ray and optical/IR emission from these systems. 
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